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Abstract 

This paper studies the correction of electromagnetic vector vertex of quarks due to the strong 
interactions. 
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Since quarks carry both electric charges and color charges, an interesting question will be 
how the strong interactions disturbs quarks' electromagnetic vector vertex. For electrons, 
the tree-level electromagnetic vector function is a simple 7^, thus the corresponding "Dirac" 
magnetic moment, which is usually expressed in term of the Landc g-factor, is: 2. However, 
higher order Feynman diagrams will differ the e-m vector vertex from a pure 7^, then g will 
slightly differ from 2. This is well known as the "anomalous magnetic moment" of electron. 
Composite particles often have a much huger anomalous magnetic moment such as proton. 
This indicates a great difference of e-m vertex function from simple 7^. 

In ordinary state free quarks do not exist, however it may exist in strange matter under 
special environments such as in neutron stars. So, it will useful to consider the "anomalous 
magnetic moment" of quarks. This is essentially to calculate the c-m vertex function of 
quarks with both higher order e-m feynman diagrams and effects of the strong interactions. 
The higher order e-m feynman diagrams are always considered as small perturbation, so we 
can safely concentrate to latter contribution. In this paper I calculate the e-m vertex func- 
tion of free quark with a relatively simple model "GCM" to describe the strong interactions, 
and employ Dyson-schwinger equation as a non-perturbation method for strong interactions 
calculation, (the Euclidean space field formulation is also employed for calculation con- 
venience in future studies; The question of how one may proceed between Euclidean and 
Minkowski space is discussed by[l].) 

The GCM generating functional of massless quarks in Euclidean space can be written as 

Z[f],f], An] = J VqDq ei"^'^*'^ q{x)[^-d+^^A^(x)]q(x)-jd'^xd'^y y^{x)Du,^{x-y)jS{y)+J d^x {fjq+QV)} 

(1) 

here j'^{x) denotes the color octet vector current j^{x) — q{x)ji^-^q{x) and Di^i,{x — y) — 
5i^i,D{x — y) denotes the gluon two-point function. 
The coordinate space vector vertex is: 

P . ^ _ 5^V{A,q^,q^') 

,[x,, X2, y) - sAniy)5q^ix,)5q^ix,) A,=o 

Where r(A,q'^\q'^'') denotes the effective action. From the relation: 

5'^T{A,q^\q''^) 



Sq'^''{xi)Sq'^''{x2) 



Eq.(2) can be expressed as: 



rJxi,X2,y)^^^r^^ (3) 
UK 1, 2,yj 5Af,{y) A^=o ^ ' 
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here g ^[Afj\ denotes the inverse quark propagator under external field and can be ex- 
panded in powers of as follows: 

= G-\xi, X2) + j d^yT^{xi, X2, y)A^{y) + ... (4) 

Which leads to the formal expansion 

g[A^]^G- ld'yGA^r,G + ... (5) 

The dressed quark propagator has the decomposition in momentum space: 

G-\p)^i^-pA{p') + B{p') 

where the functions A(p^) and B{p^) are determined by the rainbow Dyson-Schwinger equa- 
tion. 

In momentum space, Eq.(4), Eq.(5) lead to: 

Go'^9o'M-rlA, 

and 

G^g[A^] + GA^r^G 
Insert them to the rainbow equation 

and up to the first order, one gets: 

r.(p, q)-K-lJ ^A^iP - khMk + q/2)r,ik, q)G{k - (6) 
should be constrained by the vector Ward-Takahashi identity: 

q,r,{p,q) = G-\P+^-)-G-\p-^-) (7) 
Since transforms as a vector, and contains parameters p and q, it can be written as[3]: 
r^b, q) = Ai(p, q)-fi, + A2(p, q)Pi, + Aab, q)qn (8) 



Insert Eq.(8) to Eq.(7): 



Take traces of Eq.(9) of both sides: 



A2P • g + Asg' = + - - |)^) (10) 



2/ y 2 

Dot on both sides of Eq.(6) one gets: 



Dot 7^ leads to: 



4 f d^k 

= Pm7m - 3PM y (^^^(^ ~ k)-i^GiVf,G2lu (11) 



7^r^ = A17J + A27;,p^ + A37^g^ 

^ " / (0^^^ " ^)7.Gir,G'27. (12) 
To express above equations clearer I first evaluate expression 

7.G(^+|)r^(^,g)G(^-|)7. 
put expression of G{p) in it and do some calculations: 

1uG{pi)T^,{k,q)G{p2)'^y 

= 7i/^ ^/ 2N I T3( 2n [^17m + A2A:^ + Agg^]^ ^/ 2N i 2N 

«7 • Pi^(Pi) + B{pi) Z7 • P2^(P2) + B{pi) 

^ {/ + // + ///} (13) 



\p\A^{pI) + 52(p2)][p2^2(^2) + B^ipl)] 

Here pi = A; + |,P2 = ^ — f 



I = Ai{2A(p?)A(p^)(^ .p2)^^(^ .pi) - 4iA(p?)S(p^)pi^ 

- AtA{pl)B{pl)p2, - 2B{pI)B{pI)^,} (14) 



// = K2k^{-AA{pl)A{pl)p, ■ p2 + 2iA{pl)B{pl) (7 • p,) 

+ 2iA(pl)B{pl)(^.p2)+mPl)B{pl))} (15) 
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Ill = K^q^{-AAipl)A{pl)p, ■ p2 + 2iA{pl)B{pl){^ ■ p,)} 

+ 2tA{pl)B{pl){^-p2)+AB{pl)B{pl))} (16) 



In order to have a solution of function, insert expression (13) to Eq.(ll), Eq.(12) and 
take traces on both sides: 



VA2 + 4p-gA3 = --p,Jj^^D{p-k)tr{Eq.{13)} 



I 



d'^k 
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(27r)4^^^ [pjA^ipD + B^ipDUA^ipl) + B^ipl)] 



* { [-iA{pl)B{pl)p ■ pi - iA{pI)B{pI)p ■ p2\ Ai 
+ p-k \-A{pi)A{pl)p,-p2 + B{pI)B{pI)] A2 



+ p-q [-A{pI)A{pI)pi ■ p2 + B(pI)B{pI)] A3} 



(17) 



and 



I6A1 



4 

16- 3 



16-- 



d^k 

(2^ 
d^k 



D{p-k)tr{^^Eq.{13)} 
D{p-k)- 



(2vr)4^ \p\A\p\) + B\p\)\^lA\p'^ + B^{jy\) 

*{ - \^{A{p\)A{pl)p^ ■ p2 + 2B{pI)B{pI)]A, 

+[(8iA(pl)B(pl)k ■ pi + 8iA(pl)B(pl)k ■ P2)]A2 

+ [{8zA{pl)B{pl)k-p, + 8tA{pl)B{pl)k-p2)]A3} 



(18) 



These two equations above combined with Eq.(lO) and the rainbow Dyson-Schwing equation 
of functions A(p^) and -B(p^) 



[A{p') - ly 
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8 f d^q 
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A{q^)p ■ q 



q^A^{q^) + B^{q^) 
B{q') 



(19) 



will in principle fully restrict the nonperturbative vector vertex of quarks. By numerically 
study Eq.(lO) Eq.(17) Eq.(18) and Eq.(19) with appropriate gluon two-point function, one 
can get the nonperturbative vector vertex and quark propagater in Euclidean space. 

In order to have a qualitative understanding of this method, a particularly simple model 
gluon two-point function is used as follows: 



D{p-k)^g5^{p-k) 
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(20) 



when g <^ p, we have j9 + g/2 ^ p — g/2 ^ j9, Eq(lO), Eq(17), Eq(18) have a simple solution 
of Aj with parameters A{p'^) and B{p^) : 

_ 127r^[E^ + 37r^(AV + B^f/g - AY] 
^ ~ g + 367r8(A2p2 + B^y/g - ISvrMV + Gn^B^ 

2AiAB7r^ 

A2 - 



A, 



g + 367r8(AV + 52)2/^ - I^tt^A^ + 671^52 
2AiABTT^p ■ q 



g + 367r8(A2p2 + 52)2/^ - 15tt*AY + Qt^^B^ 



g2 



(21) 



Here A, B are solutions of the rainbow equation(19) 

Insert the same two-point function of gluon to Eq.(19): 



S(/)=0, = + jLj)"'^] (22) 



B(/) = 0. = + + (23) 



and: 



B{p') = - A{p') = 2 for p'<^, (24) 



Solution(24) generates a dynamical quark mass and breaks chiral symmetry sponta- 
neously. Insert this solution to Eq.(21), I get the e-m vertex function of quarks: 

r - - ^ 16i 16ip ■ q 

7^^ + - Ap^y/^P' + - 4^2)1/2^2 l^^J 
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